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SUMMARY 

A computer program Is developed for calculating the statistical fast frac- 
ture reliability and failure probability of ceramic components. The program 
Includes the two-parameter Welbull material fracture strength distribution 
model, using the principle of Independent action for polyaxlal stress states 
and Batdorf's shear-sensitive as well as shear-insensitive crack theories, all 
for volume distributed flaws In macroscoplcal ly Isotropic solids. Both penny- 
shaped cracks and Griffith cracks are Included In the Batdorf shear-sensitive 
crack response calculations, using Griffith's maximum tensile stress or crit- 
ical coplanar strain energy release rate criteria to predict mixed mode frac- 
ture. Welbull material parameters can also be ca^ulated from modulus of 
rupture bar tests, using the least squares method with known specimen geometry 
and fracture data. The reliability prediction analysis uses MSC/NASTRAN 
stress, temperature and volume output, obtained from the use of three- 
dimensional, quadratic. Isoparametric, or axlsymmetrlc finite elements. The 
statistical fast fracture theories employed, along with selected Input and 
output formats and options, are summarised. An example problem to demonstrate 
various features of the program Is Included. 


INTRODUCTION 

The attractive physical and mechanical properties of modern ceramics - 
high temperature strength, light weight, excellent erosion, corrosion and oxi- 
dation resistance, low thermal conductivity, low cost, and wide availability - 
have made ceramics an Increasingly Important structural material. The poten- 
tial of ceramics In demanding structural applications Is especially attractive 
when resistance to hlgn temperatures, such as In heat engines, Is the main 
concern. With today's needs for more fuel efficient transportation, multifuel 
engine capability and reduced emissions, advanced engines, operating at much 
higher temperatures with ceramic components, appear to be mandatory to econom- 
ically meet these national objectives. 

However, ceramics, like all other brittle materials, display linear 
stress-strain behavior from zero to fracture. The lack of ductility and yield- 
ing capability give ceramic materials their most undesirable characteristics 
such as low strain tolerance, low fracture toughness, and large variation In 
observed fracture strength. This wide variation of material strength Is due to 
the nature and distribution of Intrinsic microscopic flaws, which are unavoid- 
ably present as a result of materials processing operations. Failure In 
ceramics usually Initiates at a single weakest flaw when the local stress there 
reaches a critical value. 8ecause of the large scatter In strength, designers 
today use statistics and reliability analysis for the failure prediction of 
brittle material components, which may be suDject to arbitrary loadings and 
multidimensional stress states. 


The first probabilistic approach used to account for the seal r In frac- 
ture strength of brittle materials was Introduced by Welbull (ref. .). His 
analysis was based on the weakest link theory (WLT) and assumed a unique 
strength cumulative distribution for the uniaxial fracture data obtained from 
simple specimen tests. Welbull also proposed a method for calculating the 
failure probability In multidimensional stress fields when using material 
parameters obtained from uniaxial tests. His approach basically Invol es cal- 
culating the risk of rupture by averaging the tensile stress In all directions. 
This Is Intuitively plausible but not rigorous, and consequently other models 
were Introduced. The most widely used among them Is the assumption that prin- 
cipal stresses act Independently (PIA) (ref. 2). This Is a very convenient 
formulation because of Its simplicity. However, as shown by several Investi- 
gators these models can lead to unsafe estimates of failure probability, since 
they both ?glect the shear force, and In case of the PIA hypothesis, the 
effects of combined local principal stresses (ref. 3). 

An Important element of failure predictive theories Is the crack extension 
criterion. In the classical Welbull formulation, a normal stress criterion Is 
used, which Is likely to be correct when the dominant crack Is normal to a 
uniaxial tensile stress. However, In a multlaxial stress field with flaws 
orientated at arbitrary angles to the applied stresses both normal tensile 
stresses and In-plane shear stresses will Influence the deformation and frac- 
ture processes (refs. 4 to 6), and lead to a different fracture response than 
that of the uniaxial case. Several fracture criteria have been proposed 
(refs. 7 to 9), with the critical coplanar strain enercy release rate, G c , 
criterion, among those available In this study, leading to the best agreement 
with available brittle material experimental data. 

The primary objective cf this report Is to develop a public domain com- 
puter program which will be coupled to a general purpose finite element code, 
such as MSC/NASTRAN (ref. 10), to predict the fast fracture failure probability 
of ceramic components due to the presence of volume type flaws. The user Is 
given various options to select currently available fracture models In addltlci 
to calculating statistical material parameters. Two versions of the program 
are presently available which are designated as SCARE1 (Structural Ceramics 
Analysis and Reliability Evaluation) and SCARE2, respectively (ref. 11). In 
SCARE1 , the finite element centroldal principal stresses are taken as constant 
throughout the element volume and the convergence of the mesh for accurate 
stress analysis leads to convergence for volume type flaw reliability analysis. 
However, previous results from higher order Isoparametric finite elements with 
permissible Internal stress gradient 1 " showed that the accuracy of failure pre- 
dictions Is significantly Improved when the finite element volumes are further 
subdivided. In the SCARE2 version, all 6-slded HEXA MSC/NASTRAN elements are 
further discretized Into 27 subelements, which are then used with Interpolated 
principal stresses to perform all analysis. 


PROGRAM CAPABILITY AND DESCRIPTION 

The basic computational elements of the post-processor program for the 
reliability analysis of structural ceramic components using the SCARE2 version 
of the code are shown In figure 1. Clearly the NASTRAN part Is totally Inde- 
pendent of SCARE, and output data from other general purpose analysis programs 
could also be used as long as similar elements are available for the thermal 
and stress analysis of the structure. For computational efficiency, all the 
experimental fracture stresses. If used to calculate material parameters, as 
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well as the elemental principal stresses are normalized. Initially, the nor- 
malizing stress used Is the average of all experimental fracture data when 
specimen test results are available. When material parameters are known, an 
appropriate value of the Welbull scale parameter, o 0 , Is used to normalize 
all stresses. 

As shown In figure 1, executing the SCARE program requires FORTRAN logical 
units 1, 3, 4, 5, and 7 In addition to those generally used In performing 
NASTRAN analysis (such as units 5, 6, and 14). The program uses these tape 
drives for Intermediate storage of large amounts of data, so that the transfer 
of Information from NASTRAN to SCARE Is done Internally rather than manually. 
MSC/NASTRAN Bulk Data Including required nodal temperatures, and element 
centroldal as well as selected corner node stresses, are stored on logical 
unit 7 where NASTRAN punch files are saved for access and further processing In 
SCARE. Since element volumes are obtained through a NASTRAN Parameter call, 
the required volume data Is taken from the printout files on unit 6 and stored 
on logical unit 3. The Input to SCARE Is handled through logical unit 5, but 
the output from SCARE had to be plated on unit 1 to avoid potential terminal 
problems. In addition to NASTRAN output files and analysis data. Input to 
SCARE Includes control Indices specifying various fracture models, temperature 
dependent material parameters If available, specimen geometry and ordered (In 
ascending order for a given temperature) fracture strength data when required 
statistical parameters are Internally calculated. 

In order to use WLT, no principal compressive stress Is permitted to 
exceed three times the maximum principal tensile stress In absolute value. If 
this criterion Is violated In an element, compressive stress state predominates 
and the corresponding reliability Is set equal to unity. Additionally, when 
using the PIA model In conjunction with Welbull stalstlcs, only tensile prin- 
cipal stresses can contribute to failure and fracture due to compression Is 
Inadmissible. 

The program has broad capabilities by allowing the user to specify temper- 
ature dependent statistical material parameters, several crack configurations 
and four fracture criteria. Uniaxial fracture data along with specimen geom- 
etry from four point modulus of rupture (MOR) bar tests can be used to cal- 
culate Welbull parameters and the Batdorf crack density coefficient. Figure 2 
contains the flowchart summarizing the available options Is fracture criteria 
and flaw configurations used to model the volume Imperfections. Note that two 
of the failure criteria are for shear-insensitive cracks, even In polyaxial 
stress states. The other two criteria are used for the more general shear- 
sensitive model. The available three crack configurations Include the spher- 
ical void, which Is Isotropic or direction Independent, and Is Inherently 
assumed In the Welbull PIA and normal stress failure theories. However, Imper- 
fections In high density, sintered ceramics are best represented by the struc- 
tural response of penny-shaped and Griffith cracks. Among the available 
criteria and crack configurations shown. In figure 2, the penny-shaped crack 
(PSC) with the G c criterion gives the highest failure probability for a 
given case, while the PIA approach yields consistently the lowest failure 
estimate. It should also be noted that the Batdorf shear-insensitive fracture 
model, although Iri slightly different form, gives Identical results to the 
originally proposed Welbull normal stress averaging method. 
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INPUT INFORMATION 


NASTRAN (NASA Structural Analysis) Is a large, comprehensive, general 
purpose finite element computer code for structural analysis, which was devel- 
oped under NASA sponsorship to fill the need for a universally available 
analysis program. In addition to the government supported version, there are 
several, greatly enhanced, proprietary versions of this program, the most 
widely known of which Is called MSC/NASTRAN. This program Is used throughout 
the world In large corporations, government laboratories, and most commercial 
data centers. The SCARE program utilizes results from only a very small frac- 
tion of available NASTRAN analysis capability. Since fast fraccure mechanical 
design of ceramic components requires only the temperature and stress distribu- 
tions, static analysis results from rigid formats 61 and 47 (In case of cyclic 
symmetry) are most often used. Ceramics are also extremely sensitive to geom- 
etric discontinuities, requiring the use of Isoparametric three-dimensional 
and quadratic axlsymmetrlc finite elements. Within MSC/NASTRAN, these elements 
are denoted as HEXA, PENTA, and TRIAX6, respectively. Although the mldedoe 
nodes of HEXA elements In NASTRAN are optional, their use when analyzlr -.ch 
SCARE Is required. It Is assumed here that analysts using the SCARE pr 3 .am 
would be fully familiar with MSC/NASTRAN, and Its Input requirements In creat- 
ing the Executive Control, Case Control, and Bulk Data decks. Preceding the 
entire NASTRAN Input file Is the system operating JCL (Job control language) 
set of commands, which usually Identify the Job, user, set time and memory 
requirements, and define NASTRAN Input, punch, plot, and printout files. If 
the self-contained, solid modeling processor, called MSGMESH, Is used to dis- 
cretize the structure, the SCARE program Includes a number of sorting routines 
to permit arbitrary numbering of elements and nodes. Figure 3 shows the 
arrangement of a typical NASTRAN Input file. 

The MSC/NASTRAN program at the Lewis laboratory runs on the CRAY 1 -S/2200 
computer In a batch-mode. Input and output to and from the CRAY Is handled 
through an IBM 370 mainframe computer, which serves as the front-end processor 
for the CRAY. Consequently, an additional set of JCL commands Is required to 
handle the Involved data sets, compiler, FORTRAN logical units, and execution 
commands. Both versions of the SCARE source program are permanently stored on 
the IBM In the user's library, where ell reliability calculations are even- 
tually made. Details of the system JCL for executing MSC/NASTRAN first, and 
then the SCARE program can be found In reference 11. These Instructions are 
unique to the computer system existing at Lewis, but are representative of the 
required commands In performing reliability analysis at other Installations. 

SCARE input requirements can be grouped Into three categories. The first 
category, called Master Control Deck, defines control Indices, Information on 
the finite element mesh and some Integer data describing crack configurations, 
fracture criteria and material parameter format. Figure 4 shows the details of 
the required Information, with explanatory notes and size limitations available 
In the program user's manual (ref. 11). The second category, called the Speci- 
men Deck, summarizes fracture specimen results needed In calculating statisti- 
cal fracture parameters, or direct material properties. Including Poisson's 
ratio, when available. There are five different entries required In this 
category. The first entry Includes the material Poisson's ratio, which Is used 
In the reliability calculations when PSC's are selected for volume Imperfec- 
tions. The second entry defines the MOR specimen geometry, which was used In 
generating fracture data. The third entry Includes experimental, extreme fiber 
fracture stresses, arranged In ascending order for a given temperature. For 
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multiple temperature tests, the temperatures must also be ordered according to 
ascending values, since calculated material parameters are Interpolated within 
SCARE. At a specified temperature, fracture readings must be unique and mul- 
tiple values of Identical magnitudes are not permitted. The number of avail- 
able fracture readings for all temperature tests must be the same. The fourth 
entry Is used when material statistical parameters are directly available as a 
function of temperature. The three required parameters are the Welbull modulus 
or shape parameter, m, the Batdorf normalized crack density coefficient or flaw 
parameter, kg, and the Welbull scale parameter, o 0 . These material param- 
eters must be so arranged that t.iey correspond to ascending order of discrete 
temperatures . The last entry In the Specimen Deck category lists the discrete 
temperatures at which fracture data or material parameters are known. Addi- 
tional explanation of the required Input, Including size limitations, can be 
found In reference 11. 

The last SCARE Input category, called the Structures Deck, contains 
results of the finite element structural analysis needed for failure probabil- 
ity predictions. These Include element volume, element or nodal temperatures, 
and element principal stresses along with the appropriate Identification 
numbers. In the present version of SCARE, which relies on MSC/NASTRAN output 
files, all of this data Is internally manipulated through subroutine ELEM, and 
the Structures Deck requires r.o specific Input by the user. It Is this Input 
data, however, that has to be carefully catalogued If another general purpose 
program were to be used or element data were to be directly read. Instead of 
using temporary storage devices. 


OUTPUT INFORMATION 

The results of MSC/NASTRAN thermal analysis are the grid point tempera- 
tures, which can be obtained at transient or steady-state conditions. After 
solution of the component temperature distribution, the most severe thermal 
gradients can be selected and combined with th* mechanical loads to obtain an 
elastic solution {rigid formats 47 or 61). The usual output from these r 1g1d 
formats Includes the nodal displacements along the global coordinate directions 
and the element stresses. Depending on the element type, normal and shear 
stresses In the local element, or material coordinate system are always aval - 
able at element centroids. Additionally, element principal stresses are cal- 
culated there for the HEXA and PENTA elements. Corner node stresses are also 
printed for these elements. For available stress recovery options, users 
should consult reference 10 and the appropriate program manuals. In addition 
to the displacements and stresses, useful parameters such as element volumes, 
element areas, component center of gravity, moments of Inertia, etc. can be 
calculated through the Parameter call feature of the program. For volume flaw 
reliability analysis, the element volumes are esst tial, since In the weakest 
link model, Integration of the stress distribution over the material volume Is 
needed. 

The first part of all SCARE output data contains an echo of Important 
NASTRAN finite element analysis results. Identifying labels, element type, 
and number of elements In the model are noted. A table of element centroldal 
principal stresses with appropriate element Identification numbers Is given. 

For the SCARE2 version, the 27 centroldal subelement principal stresses within 
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each HEXA element are listed. Element volumes and calculated element tempera- 
tures are summarized In another output table. Next, the selected fracture 
model Is Identified and the room temperature (70 # F) statistical fracture 
parameters are listed. Additionally, a table of discrete temperatures with 
corresponding material parameters, which were either Internally calculated or 
directly supplied, Is printed. For the shear-sensitive fracture models, the 
crack shape Is Identified along with a more specific description of the frac- 
ture criterion. The last table In the SCARE output file contains an element 
results summary, listing the element number, corresponding element survival 
probablll les and material parameters. Finally, the overall component proba- 
bility of failure as well as the component probability of survival are printed. 


THEORY 

The statistical nature o c fracture In engineering materials can be viewed 
from two distinct and extreme .nodels. The first was presented by Welbull and 
Is termed the weak link model. With It a structure Is characterized as a 
series of links connected In such a manner that the structure falls whenever 
any of the links fractures. In contrast the second model Is referred to as 
the bundle or parallel model for which failure Is defined only when all links 
In parallel have fractured. Structural ceramics have been observed to approach 
the weakest link hypothesis and fall when the stress Intensity factor at any 
one flaw reaches a critical value. In view of Its pessimism, WLT design Is In 
most cases conservative. Other Important features of WLT are that It predicts 
size effect and that failure of a complex component may not be Initiated at the 
point of highest nominal stress. A particularly severe flaw may be located at 
somewhat less highly stressed point and may still be the first crack to become 
critical. It Is for this reason that the entire field solution of the stresses 
must be obtained and examination of the most highly stressed point, as In 
ductile materials. Is no longer adequate. 

Experimental fracture strength data obtained from uniaxial ly loaded simple 
specimens, when arranged In ascending order, can be represented In two differ- 
ent forms. The probability density function of a random variable (fracture 
strength) Is a mathematical function that best represents the data In a rela- 
tive frequency histogram, that Is failure stress versus number of failures. 

The result, typically, would be a bell shaped curve (ref. 12). Alternatively, 
fracture data can be plotted as stress versus failure probability which leads 
to an S-shaped curve, called the cumulative distribution function of this 
random variable. Various distribution curves have been used to characterize 
the material's fracture property. The two most commonly used distribution 
functions are the Gaussian (normal) distribution and Weibull's distribution. 

The Welbull distribution Is selected to characterize ceramic strength varia- 
tions, since che Gaussian distribution is Intrinsically associated with the 
bundle model and Is Incompatible with WLT. Consequently, the uniaxial fracture 
data Is approximated by the 3-parameter Welbull distribution, defined by 
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where Pf Is the probability of failure, o 0 Is the scale parameter with 
dimensions of stress x (volume)!/ 1 * 1 , o u Is the threshold stress which Is 
usually taken as zero, m Is the Welbull modulus which measures the degree of 
strength variability, a Is the applied tensile stress and V the stressed 
volume. 


In the analysis of failure of brittle materials subject to multlaxlal 
stress states, the Welbull model, when combined with the PIA hypothesis, yields 


P 


f 


. 1 



( 2 ) 


where o-| , ©2» an d °3 are the Principal stresses and It was assumed that 
o u * 0. Equation (2) has been widely used In the past to estimate failure 
probabilities of ceramic structures (ref. 13). The failure probability using 
the normal tensile stress averaging method, as proposed early by Welbull 
(ref. 1), and described later through an Integral formulation (ref. 14) can be 
calculated from 


1 - exp - f 


fwp f oJJ dAjdV 


(3) 


where k wp Is the polyaxlal Welbull crack density coefficient given by 



( 4 ) 


This constant can be obtained by making the result of Integrating equation (3), 
using the normal stress o n distribution on an arbitrary plane, obtained from 
the Cauchy Infinitesimal tetrahedron In principal stress space as shown In 
figure 5, for uniaxial stress cases, agree with the results obtained from the 
uniaxial, 2-parameter Welbull equation. The area Integration Is performed on 
the surface of the unit sphere where the normal stress Is tensile and neglect- 
ing regions where the normal stress Is compressive. The crack-like flaws can 
then be regarded as located In these arbitrary planes which are tangent to the 
sphere and are acted upon by o n which Is Induced by the principal stresses 
<* 1 , o 2 . a nd °3- Since equation (3) Is just the shear-insensitive case of the 
more general Batdorf (ref. 15) polyaxlal stress fracture model. Its SCARE 
Implementation follows a somewhat different format. The polyaxlal Welbull 
equation has also been extensively used In the past (ref. 14), but since It 
neglects the effects of shear loads. It also underestimates failure for the 
more general loading condition. 

In the previously described two multi-dimensional stress fracture models 
no direct use was made of the hypothesis that fractures are due to crack 
growth. In references 4 and 15, attention Is focused on the cracks and their 
failure under stress. Since there Is not as yet a consensus regarding how to 
treat mixed mode fracture, even In ductile materials, the SCARE program 
Includes several fracture criteria and flaw shapes. Rufln et al. (ref. 16) 
recently compared results obtained from various fracture models and experi- 
mental tests, with similar work being reported In reference 5. 
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Consider now a small uniformly stressed material element of volume av. 
The probability of failure under an applied state of stress can be written as 
(ref. 16) 


P f = Pt P 2 


(5) 


where Pi Is the probability of existence In aV of a crack having a cri- 
tical stress In the range of o cr to o cr + do cr , and P? denotes the proba- 
bility that a crack of critical stress o cr will be oriented In a direction 
such that an effective stress o e equals or exceeds o cr . u cr Is defined 
as the remote, uniaxial, normal fracture stress of a given crack. Failure will 
occur when the effective stress (a function of chosen crack configuration and 
fracture criterion) exceeds o c - for a particular crack. P-| has the form 


P 1 " * V do 


dN 


cr 


(a ) do 
' cr' cr 


( 6 ) 


and 


where N(o cr ) Is the crack density function (the density of cracks having a 
critical stress < o cr ) and fi Is the solid angle In principal stress space 
containing all the crack orientations for which o e > o cr . Using the weakest 
link theory, the overall failure probability can be calculated from (ref. 15) 



The crack density function N(o cr ) Is a material constant and Is Inde- 
pendent of stress state. It Is usually expressed as a power function of o cr , 
that Is N(o cr ) = kgof r , where the flaw distribution parameters kg and m 
can be evaluated from experimental data using uniaxial or equlblaxlal tension 
specimens. Batdorf (ref. 15) Initially proposed a Taylor series expansion for 
K(o cr ), but this method had computational difficulties. Recently, a more 
convenient Integral equation approach was formulated and extended to the use 
of data from 4-polnt MOR bar tests (ref. 16). 

The statistical analysis of fracture Is greatly simplified by assuming 
that cracks are shear-insensitive. For this case fracture occurs when 
<j e = o cr = <* n and there Is no need to specify the crack shape or the mate- 
rial's Poisson's ratio. Note that the crack size Is never used In statistical 
fracture theories and Is always eliminated from the analysis. Since for uni- 
axial loading shear-insensitive cracks are assumed to dominate the fracture 
process, a e Is defined such that In the absence of shear on the crack plane, 
a e = o n . In a similar manner, when the G c criterion Is used, we define 
o e as the uniaxial normal stress that would Induce the same energy release 
rate as the actual stress. The same Ideas can be extended to noncoplanar crack 
growth criteria, such as the maximum G or strain energy density, to define 
o e for those applications. In any event, for polyaxlal stress states, the 
effective stress o e Is a function of both o n and r, where t Is the shear 
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stress In the crack plane. Similarly to calculating o n at any point, the 
Cauchy Infinitesimal tetrahedron of Fig. 5 can also be used to obtain t on 
the same arbitrary crack plane. Reference 4 gives effective stress expressions 
for two crack shapes using the maximum tensile stress and G c fracture cri- 
teria. The same four options for shear-sensitive cracks are available In 
SCARE, a- can be seen In figure 2. The best choice among them Is o e for 
PSC's, given by 


°e = t2/(1 ‘ °- 5 w)2 (9) 

where v Is Poisson's ratio. 

The solid angle Si depends on the fracture criterion selected, the 
assumed crack configuration and on the applied stress state. Closed form 
expressions for Q can be derived for analytically simple fracture criteria 
In uniaxial and balanced biaxial stress states (ref. 15). Assuming a uniaxial 
stress, o, and the normal stress (shear-insensitive) fracture criterion, we 
obtain (ref. 5) 


P 




( 10 ) 


Note that when a shear-sensitive fracture criterion Is used, the crack shape 
must also be specified. In general for three-dimensional stress states. Si 
must be determined numerically. Using the shear-insensitive case as an 
example, we obtain at fracture (ref. 15) 


2 2 2 

o = o = o = o- + (o, - a. ) COS a + (o„ - a„) COS 0 sin a (11) 

6 n c r j i o c j 

where from figure 5 angles a and 0 define the crack plane In principal 
stress space, on which o n and t act. Using direction cosines 1, m, 
and n, the equilibrium of forces on the Cauchy tetrahedron yields values of 
o n and t In terms of the principal stresses and the angles a and 3 
(ref. 11). It Is computationally convenient to define <p = cos? 0. Then 
equation (11) can tic rewritten as 


a.|<p + a^tp + a^ = 0 


( 12 ) 


where a-| = 0, a? = (02 - 03 ) sin? a and a_ = (o-| - 03 ) 
cos? a + 03 - o cr . Solving for <p gives 

“ a 3 °cr " °3 " (o l ' °3 ) C0$2 a 
9 = — = - 
2 (o ? - o 3 ) s 1 n‘ a 

If we define 0 = cos -1 V<p, thtr. °2 can be calculated from (ref. 15) 


(13) 
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suits m equation (1) are compared to the expression given . nuatlon 
(15), ano the relationship between o 0 and C Is derived by making Pf from 
the two equations agree. For a rectangular beam of width w, length L-| 
between symmetrically placed outer loads and length l? between Inner loads, 
the result Is 


a 


o 



(L 1 ♦ ml_ 2 ) 
C(m ♦ l) 2 


1/m 


(16) 


In addition to obtaining o 0 and m, the SCARE program requires knowledge 
of kg. We can evaluate kg from 4-polnt flexure data by substituting 
equation (10) Into equation (8). Integration by parts of the results gives 


P $ » exp 



"<°cr> 



dV 


(17) 


We again utilize the power function form of N(o cr ) In equation (17) and 
carry out the stress Integration. Similarly to the Welbull analysis, o-| Is 
expressed In terms of of and the beam height. The volume Integration Is 
then performed over the tensile portion of the beam. Including effects of 
changing o-j along the beam length. Results from this Integration are com- 
pared to equation (15) and the relationship between kg and C Is derived by 
setting Pf equal from the two equations. For a 4-polnt loaded beam speci- 
men having a rectangular cross section, we obtain 

k B ■ < 2ra * ] (1 


By comparing equations (16) and (18), we conclude that when using the normal 
stress failure criterion, kg and a 0 are related by 


k 0 = (2m ♦ 1) 



(19) 


However, equation (19) changes when the G c criterion Is used (ref. 5). 


EXAMPLE 

In order to validate SCARE, a number of example problems were analyzed 
from the open literature (ref. 11). Among them failure probability prediction, 
were made, using Batdorf's shear-insensitive fracture model, for a silicon 
nitride disk rotating at various angular velocities (ref. 14). The dimensions 
of the disk along with appropriate material statistical parameter data are 
given In figure 6. Because of the simple geometry, only eight HEXA elements 
were used In one 15° sector MSC/NASTRAN model of the disk. The calculated 
NASTRAN stresses and volumes both were within approximately 1 percent of the 
available closed form answers. Both SCARE1 and SCARE2 predictions were gen- 
erated and the results were compared to those listed In reference 14. Reli- 
ability calculations were also made at various speeds using other fracture 
theories. Selected results from these ^ijialyses are shown In figure 7 and 
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In addition to obtaining o 0 and m, the SCARE program requires knowledge 
of kg. We can evaluate kg from 4-polnt flexure data by substituting 
equation (10) Into equation (8). Integration by parts of the results gives 


P $ * exp 





dd 


cr 
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We again utilize the power function form of N(o cr ) In equation (17) and 
carry out the stress Integration. Similarly to the Welbull analysis, o-| Is 
expressed In terms of of and the beam height. The volume Integration Is 
then performed over the tensile portion of the beam. Including effects of 
changing o-| along the beam length. Results from this Integration are com- 
pared to equation (15) and the relationship between kg and C Is derived by 
setting Pf equal from the two equations. For a 4-polnt loaded beam speci- 
men having a rectangular cross section, we obtain 


k 0 = (2m ♦ 1) 


2C(m + l) 2 
wh(L 1 + ml_ 2 ) 


(18) 


By comparing equations (16) and (18), we conclude that when using the normal 
stress failure criterion, kg and o 0 are related by 


kg = (2m + 1) 



(19) 


However, equation (19) changes when the G c criterion Is used (ref. 5). 


EXAMPLE 

In order to validate SCARE, a number of example problems were analyzed 
from the open literature (ref. 11). Among them failure probability prediction, 
were made, using Batdorf's shear-insensitive fracture model, for a silicon 
nitride disk rotating at various angular velocities (ref. 14). The dimensions 
of the disk along with appropriate material statistical parameter data are 
given In figure 6. Because of the simple geometry, only eight HEXA elements 
were used In one 15° sector MSC/NASTRAN model of the disk. The calculated 
NASTRAN stresses and volumes both were within approximately 1 percent of the 
available closed form answers. Both SCARE1 and SCARE2 predictions were gen- 
erated and the results were compared to those listed In reference 14. Reli- 
ability calculations were also made at various speeds using other fracture 
theories. Selected results from these ^alyses are shown In figure 7 and 


table II. For a given speed, the Welbull PIA model yields clearly the lowest 
failure estimate, while the shear-sensitive PSC w'th the G c criterion gives 
the highest. The agreement between SCARE2 results and those In reference 14 
was within 10 percent, with the difference probably due to the different 
stress-volume data used In solving the reliability problem. Laboratory mea- 
surements agree best with the selected shear-sensitive fracture model, as car. 
be noted, especially In the high failure probabilities range. Since only seven 
disks were fracture tested compared to 85 MOR specimens, there Is some concern 
about the accuracy of the experimental disk Welbull modulus of 4.95, which 
causes the greater difference between experimental and predicted Pf at 
lower failure probabilities. 


CONCLUSIONS 

A general purpose, statistical, fast fracture failure probability code 
has been generated, which Is coupled with MSC/NASTRAN, and can be used to 
design structural ceramics components. The program Includes a number of widely 
used polyaxlal fracture models, appropriate extreme value statistics and the 
ability to calculate material failure distribution parameters, all for volume 
distributed flaws. Current work Includes extension of this same capability to 
blmodal flaw populations, where failure due to extrinsic defects Is a con- 
current possibility. The addition of more advanced failure criteria which 
permit out-of -plane crack extension Is also planned. Finally, the problem of 
a transversely loaded circular plate will be Investigated, both analytically 
and experimentally, to resolve some of the contradictory trends reported In 
references 5 and 16. 
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TA8LF I. - FORMS OF P ? FOR VARIOUS SHEAR-SENSITIVE FRACTURE CRITERIA AND 
SELECTED CRACK CONFIGURATIONS (o ? * 03) 


Fracture 

criterion 

Crack 

configuration 

P 2 

Maximum tensile 
stress 

Griffith crack 
(G.C.) 

a 1 r 1 _ l 

~~x~ * - J cos V* sin 0 da 
0. * 0 



where 

s 

1 


o -a. + Ja? - 4 a. a. 

„ 2 - 1 2 1 J 

* cos 8 27j 

; 


and 

i 

i 


0 4 ! 

a l = (°2 " sin a ’ 

i 

; 


2 2 2 1 
= ( c 2 ~ sin a 1 2 ( 0 j cos a + sin a) 

; 


- 4 ' 1 cr - °3 - °?> 



2 2 2 2 
= (oj cos a ♦ J3 sin a) - 4° cr (02 COS a j 


0 2 2 2 2 2 i 

♦ 0 ^ Sin^ a) - ( 0 ^ COS a "*■ S l n a) + ^° cr | 

J 1 I j 

i 

Penny-shaped 
crack (PSC) 

2 4 * 

- o 3 ) sin a < 

! 

1 i 

i 1 

; 1 

2 2 2 ' 
&2 = Dj(a-, - 03 ) s ’n a ( 2 (oj cos a + 03 sir. a) j 

1 

4 j 

- 07 °cr - °3 * °2 ) 

1 | 

j 

i ; 

2 2 2 2 1 
a^ - D 0 *2 cos a + 03 sin a) - 4a cr (°j cos a I 

| j 

l i 

1 1 

2 2 2 2 2 2 i 

♦ 0 ^ Sin a) - COS^ a + Sin* a) + 4o cr j 

jStrain energy 
jrelease rate 

i Griffith crack 
i (G.C.) 

j 

2 2 2 2.2 
a^ r - Oj COS a - S 1 n a 

■> - 1 — 2 — ~i 

(o^ - 03 ) 5in a \ 

1 

j 

f ' — 

j Penny-shaped 
; crack (PSC) 

2.4 1 

1 a l ~ ^2 ( °2 ~ e 3 ' sin ct i 

i 

1 


a- = Dj{<;? “ 4) ' ln? a + 2D 2^ °2 _ 03) sin a a 



! ! 

0 2 

• ( a ^ sin*" a ♦ c A cos a) ! 

i 1 


2 2 2 2 2 
a 3 = Dj ( a 1 cos 0+03 sin o) + cos a 


i . ? ,2 2 

I * °3 Sln a) - °cr 

j where 

! r, 1 a -v(l - 0.25 v) 

D t = j , Up = 7- 

I (1 - 0.5 vT 1 (1 - 0.5 vT 


TABLE II. - EXAMPLE 1 FAILURE PROBABILITIES AS A FUNCTION OF ROTATIONAL SPEED FOR VARIOUS 

FRACTURE MODELS 


[m = 7.65; o 0 = 74.82 MPa (m) 0 - 39 ^; k B „S = 16.30; NGP = 5] 


Angular 

speed 

(rpm) 

SCARE 2 

(HEXA elements) 

1 SCARtl 

| (HEXA 

j elements) 

Ford (Ref. 14) 
( ax i symmetric 
elements) 

Ford (Ref. 14) 


Weibul 1 

Batdorf shear- 

: Batdorf shear- 

t" 

j Batdorf shear- 

Shear-insensitive 

Experimental 


PI A 

insensitive 

1 sensitive PSC 

1 insensitive 






! 6 r criterion 

1 

1 



70 000 

0.0021 

0.0026 

i 0.0078 

0.0022 

0.0023 

0.0583 

75 000 

.0061 

.0075 

! .0222 

1 .0064 

.0067 

.1121 

80 000 

.0163 

.0201 

! .0584 


.0170 

.0179 

.2017 

85 000 

.0412 

.0505 

1 .1426 


.0426 1 

.0446 

.3367 

93 000 

.1530 

.1850 

.4549 


.1579 i 

.1650 

.6321 

100 000 

.3954 

.4623 

.8410 


.4074 ' 

.4223 

.8714 

104 000 

.6000 

.6763 

.9649 

.6124 

.6321 

.9514 

110 000 

.8847 

.9301 

i .9996 

.8931 

.9055 

.9949 

114 000 

.9736 

.9900 

j 1.0000 

.9792 

.9830 

.9994 


TA3LE III. - EXAMPLE 1 FAILURE PROBABILITIES AS A FUNCTION OF APPLIED PRESSURE 


FOR VARIOUS FRACTURE MODELS 


[m = 28.53; t> 0 = 36 200 psi (in) 0 - 105 ; k B c.£ = 58.06; NGP = 15] 


Pressure, j 
MPa { Dsi ) 1 

i 

SCARE 2 
( three-dimensional 

! 

elements) 

i 

SCARE i 
( three-dimensional 

elements) 

j SCARE 1 

! ( ax i symmetric 
elements) 

i Weibul 1 
! PIA 

i 

.. L 

Batdorf shear- 
insensitive 

1 Batdorf shear- 
sensitive PSC 
■! en . re . ra . cr . 

Weibul 1 
PIA 

j Batdorf shear- 
! insensitive 

i 

j 

Batdorf shear- 
sensitive PSC 
en. re. ra. cr. 

; Batdorf shear- 
i sensitive PSC 
; en . re. ra. cr. 

1.31 (190)| 


0.0039 

0.0065 


! 0.0030 

0.0049 

! 

0.0064 

1.38 (200) ! 

0.0037 

.0168 

.0279 

0.0029 

1 .0128 

.0211 

.0275 

1.45 (210) 

.0149 

.0658 

.1078 | 

.0115 

1 .0504 

.0821 

i .1061 

1.52 (220)| 

.0552 

.2265 

.3495 ! 

.0425 

.1770 

.2760 

1 .3448 

1.59 (230), 

. 1828 

.5987 

. 7832 ! 

.1430 

! .4996 

.6827 

.7775 

1.66 (240) j 

.4934 

.9538 

.9942 ' 

.4054 

i .9028 

.9790 

.9937 

1.72 ( 250 ) j 

.8869 

.9999 

1.0000 i 

.8110 

.9994 

1.0000 

1.0000 

1.79 ( 260) | 

.9987 

1.0000 

! 1.0000 j .9939 

1 1.0000 

1.0000 

, 1.0000 



Fiq. l Computational elements of the SCARF2 reliability analysis program. 













MAIN PROGRAM G. C. - GRIFFITH CRACK 
"~“~I PSC - PENNY-SHAPED CRACK 



subroutine ANGLE. 





COLUMNS 

VARIABLE 

ENTRY (FORMAT - 1115) 

1-5 

1 

1 

! 

ID1 

CONTROL INDEX FOR EXPERIMENTAL DATA 

1: PURE BENDING TEST DATA 
2: 4-pt BENDING TEST DATA 
3: ALL THE MATERIAL PARAMETERS VMT, VKT 
AND VSPT ARE KNOWN AS INPUT 

6-10 

ID2 

i 

i 

CONTROL INDEX FOR FRACTURE CRITERIA 

1: SHEAR-INSENSITIVE, NORMAL STRESS 
CRITERION 

2: MAXIMUM TENSILE STRESS CRITERION 
3: ENERGY RELEASE RATE CRITERION 
4: WEIBULL PIA SHEAR-INSENSITIVE MODEL 

11-15 

103 

CONTROL INDEX FOR CRACK SHAPES 

1: GRIFFITH TYPE CRACK 
2: PENN Y-SHAPED CRACK 

16-20 

NE 

TOTAL NUMBER OF ELEMENTS IN MODEL 

21-25 

NH 

NUMBER OF HEXA ELEMENTS IN MODEL 

26-30 

NP 

| 

NUMBER OF PENTA ELEMENTS IN MODEL 

31-35 

NA 

NUMBER OF TRIAX6 ELEMENTS IN MODEL 

36-40 

NT 

TOTAL NUMBER OF SPECIMENS IN EACH SET AT 
A GIVEN TEMPERATURE 

41-45 

NGP 

NUMBER OF GAUSSIAN QUADRATURE POINTS 

46-50 

NS 

NUMBER OF SEGMENTS IN CYCLIC SYMMETRY 
PROBLEMS 

51-55 

JT 

NUMBER OF TEST TEMPERATURES AT WHICH 
MATERIAL DATA IS SPECIFIED 


Fig. 4 SCARE master control deck data requirements. 






Fig. 5 Stresses on Cauchy infinitesimal tetrahedron in principal 
stress space. 



DATA: 

NC- 132 HOT PRESSED Si^ 4 
m = 7. 65 

c 0 = 74. 82 MPa nv 3922 
k B - 16. 30 

r j = 6. 35 mm (.25 in) 
r 0 = 41.275 mm (1.625 in) 
t ■ 3.80 mm (.15 in) 

RPM RANGE - 70K TO 114K 


Fig. 6 Example 1 - rotating annular disk (Ref. 14). 



PROBABILITY OF FAILURE 



Fig. 7 Example 1 probability of failure vs disk rotational 
speed for various fracture models (SCARE2 data). 
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